Search for: All records

We conduct a systematic search for quasars with periodic variations from the archival photometric data of the Zwicky Transient Facility by cross-matching with the quasar catalogues of the Sloan Digital Sky Survey and Véron-Cetty and Véron. We first select out 184 primitive periodic candidates using the generalized Lomb–Scargle periodogram and autocorrelation function and then estimate their statistical significance of periodicity based on two red-noise models, i.e. damped random walk (DRW) and single power-law (SPL) models. As such, we finally identify 106 (DRW) and 86 (SPL) candidates with the most significant periodic variations out of 143 700 quasars. We further compare DRW and SPL models using Bayes factors, which indicate a relative preference of the SPL model for our primitive sample. We thus adopt the candidates identified with SPL as the final sample and summarize its basic properties. We extend the light curves of the selected candidates by supplying other archival survey data to verify their periodicity. However, only three candidates (with 6–8 cycles of periods) meet the selection criteria. This result clearly implies that, instead of being strictly periodic, the variability must be quasi-periodic or caused by stochastic red-noise. This exerts a challenge to the existing search approaches and calls for developing new effective methods.

 

We performed a rigorous reverberation-mapping analysis of the broad-line region (BLR) in a highly accreting (L/L_Edd= 0.74–3.4) active galactic nucleus, Markarian 142 (Mrk 142), for the first time using concurrent observations of the inner accretion disk and the BLR to determine a time lag for the Hβλ4861 emission relative to the ultraviolet (UV) continuum variations. We used continuum data taken with the Niel Gehrels Swift Observatory in theUVW2 band, and the Las Cumbres Observatory, Dan Zowada Memorial Observatory, and Liverpool Telescope in thegband, as part of the broader Mrk 142 multiwavelength monitoring campaign in 2019. We obtained new spectroscopic observations covering the Hβbroad emission line in the optical from the Gemini North Telescope and the Lijiang 2.4 m Telescope for a total of 102 epochs (over a period of 8 months) contemporaneous to the continuum data. Our primary result states a UV-to-Hβtime lag of${8.68}_{-0.72}^{+0.75}$days in Mrk 142 obtained from light-curve analysis with a Python-based running optimal average algorithm. We placed our new measurements for Mrk 142 on the optical and UV radius–luminosity relations for NGC 5548 to understand the nature of the continuum driver. The positions of Mrk 142 on the scaling relations suggest that UV is closer to the “true” driving continuum than the optical. Furthermore, we obtain$\log (M_{•}/M_{\odot })$= 6.32 ± 0.29 assuming UV as the primary driving continuum.

 

We compile optical and mid-infrared light curves from the time-domain surveys (i.e. CRTS, PTF, ZTF, and ASAS-SN) and Wide-field Infrared Survey Explorer (WISE) archive for a selected sample of active galactic nuclei (AGNs) with Hβ reverberation mapping (RM) measurements. We measure the time lags (and thus torus sizes) of W1 (∼3.4 $\mu$ m) and W2 (∼4.6 $\mu$ m) band light curves relative to the optical one using the MICA method. Through Hβ RM, the sample has well-measured AGN properties, therefore allowing us to reliably constrain the relations between torus sizes and AGN properties. We perform linear regressions for the relations between torus sizes and 5100 Å luminosities ($R\propto L_{5100}^{\beta }$ ) in two cases: β = 0.5 and β set free. The latter case yields β ≈ 0.37 ± 0.028 for both W1 and W2 bands, shallower than the expected value of 0.5, possibly due to the dependence of torus size on accretion rate. For β = 0.5, by combining with the previous K band RM measurements, we obtain the characteristic broad-line region (BLR) and tours sizes following RBLR:RK :RW1:RW2 = 1.0:6.2:9.2:11.2. We investigate the deviations of the W1 and W2 band observed torus sizes from the corresponding best-fitting relations (with β = 0.5) and find that they both are correlated with accretion rate. As the accretion rate increases, the torus sizes tend to be shortened compared to the anticipated sizes from the best-fitting relations, similar to the behaviour found in BLRs. Such behaviours can be explained by the self-shadowing effect of slim discs. This is further supported by ratios of the W1 and W2 band torus sizes to BLR sizes, which do not show significant correlations with AGN properties.

 

The AGN STORM 2 campaign is a large, multiwavelength reverberation mapping project designed to trace out the structure of Mrk 817 from the inner accretion disk to the broad emission line region and out to the dusty torus. As part of this campaign, Swift performed daily monitoring of Mrk 817 for approximately 15 months, obtaining observations in X-rays and six UV/optical filters. The X-ray monitoring shows that Mrk 817 was in a significantly fainter state than in previous observations, with only a brief flare where it reached prior flux levels. The X-ray spectrum is heavily obscured. The UV/optical light curves show significant variability throughout the campaign and are well correlated with one another, but uncorrelated with the X-rays. Combining the Swift UV/optical light curves with Hubble Space Telescope UV continuum light curves, we measure interband continuum lags,τ(λ), that increase with increasing wavelength roughly followingτ(λ) ∝λ^4/3, the dependence expected for a geometrically thin, optically thick, centrally illuminated disk. Modeling of the light curves reveals a period at the beginning of the campaign where the response of the continuum is suppressed compared to later in the light curve—the light curves are not simple shifted and scaled versions of each other. The interval of suppressed response corresponds to a period of high UV line and X-ray absorption, and reduced emission line variability amplitudes. We suggest that this indicates a significant contribution to the continuum from the broad-line region gas that sees an absorbed ionizing continuum.

 

Determining black hole masses and accretion rates with better accuracy and precision is crucial for understanding quasars as a population. These are fundamental physical properties that underpin models of active galactic nuclei. A primary technique to measure the black hole mass employs the reverberation mapping of low-redshift quasars, which is then extended via the radius–luminosity relationship for the broad-line region to estimate masses based on single-epoch spectra. An updated radius–luminosity relationship incorporates the flux ratio of optical Fe ii to H β ($\equiv \mathcal {R}_{\rm Fe}$) to correct for a bias in which more highly accreting systems have smaller line-emitting regions than previously realized. In this work, we demonstrate and quantify the effect of using this Fe-corrected radius-luminosity relationship on mass estimation by employing archival data sets possessing rest-frame optical spectra over a wide range of redshifts. We find that failure to use an Fe-corrected radius predictor results in overestimated single-epoch black hole masses for the most highly accreting quasars. Their accretion rate measures (LBol/LEdd and $\dot{\mathscr{M}}$ ) are similarly underestimated. The strongest Fe-emitting quasars belong to two classes: high-z quasars with rest-frame optical spectra, which, given their extremely high luminosities, require high accretion rates, and their low-z analogues, which, given their low black holes masses, must have high accretion rates to meet survey flux limits. These classes have mass corrections downward of about a factor of two, on average. These results strengthen the association of the dominant Eigenvector 1 parameter $\mathcal {R}_{\rm Fe}$ with the accretion process.

 

Abstract Detecting continuous nanohertz gravitational waves (GWs) generated by individual close binaries of supermassive black holes (CB-SMBHs) is one of the primary objectives of pulsar timing arrays (PTAs). The detection sensitivity is slated to increase significantly as the number of well-timed millisecond pulsars will increase by more than an order of magnitude with the advent of next-generation radio telescopes. Currently, the Bayesian analysis pipeline using parallel tempering Markov Chain Monte Carlo has been applied in multiple studies for CB-SMBH searches, but it may be challenged by the high dimensionality of the parameter space for future large-scale PTAs. One solution is to reduce the dimensionality by maximizing or marginalizing over uninformative parameters semianalytically, but it is not clear whether this approach can be extended to more complex signal models without making overly simplified assumptions. Recently, the method of diffusive nested (DNest) sampling has shown capability in coping with high dimensionality and multimodality effectively in Bayesian analysis. In this paper, we apply DNest to search for continuous GWs in simulated pulsar timing residuals and find that it performs well in terms of accuracy, robustness, and efficiency for a PTA including  ( 10 2 ) pulsars. DNest also allows a simultaneous search of multiple sources elegantly, which demonstrates its scalability and general applicability. Our results show that it is convenient and also highly beneficial to include DNest in current toolboxes of PTA analysis. 

We report the results of long-term reverberation mapping campaigns of the nearby active galactic nuclei (AGNs) NGC 4151, spanning from 1994 to 2022, based on archived observations of the FAST Spectrograph Publicly Archived Programs and our new observations with the 2.3 m telescope at the Wyoming Infrared Observatory. We reduce and calibrate all the spectra in a consistent way, and derive light curves of the broad H β line and 5100 Å continuum. Continuum light curves are also constructed using public archival photometric data to increase sampling cadences. We subtract the host galaxy contamination using Hubble Space Telescope imaging to correct fluxes of the calibrated light curves. Utilizing the long-term archival photometric data, we complete the absolute flux-calibration of the AGN continuum. We find that the H β time delays are correlated with the 5100 Å luminosities as $\tau _{\rm H\beta }\propto L_{5100}^{0.46\pm 0.16}$. This is remarkably consistent with Bentz et al. (2013)’s global size–luminosity relationship of AGNs. Moreover, the data sets for five of the seasons allow us to obtain the velocity-resolved delays of the H β line, showing diverse structures (outflows, inflows, and discs). Combining our results with previous independent measurements, we find the measured dynamics of the H β broad-line region (BLR) are possibly related to the long-term trend of the luminosity. There is also a possible additional ∼1.86 yr time lag between the variation in BLR radius and luminosity. These results suggest that dynamical changes in the BLR may be driven by the effects of radiation pressure.

 

Abstract We present reverberation mapping measurements for the prominent ultraviolet broad emission lines of the active galactic nucleus Mrk 817 using 165 spectra obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our ultraviolet observations are accompanied by X-ray, optical, and near-infrared observations as part of the AGN Space Telescope and Optical Reverberation Mapping Program 2 (AGN STORM 2). Using the cross-correlation lag analysis method, we find significant correlated variations in the continuum and emission-line light curves. We measure rest-frame delayed responses between the far-ultraviolet continuum at 1180 Å and Ly α λ 1215 Å ( 10.4 − 1.4 + 1.6 days), N v λ 1240 Å ( 15.5 − 4.8 + 1.0 days), Si iv + ]O iv λ 1397 Å ( 8.2 − 1.4 + 1.4 days), C iv λ 1549 Å ( 11.8 − 2.8 + 3.0 days), and He ii λ 1640 Å ( 9.0 − 1.9 + 4.5 days) using segments of the emission-line profile that are unaffected by absorption and blending, which results in sampling different velocity ranges for each line. However, we find that the emission-line responses to continuum variations are more complex than a simple smoothed, shifted, and scaled version of the continuum light curve. We also measure velocity-resolved lags for the Ly α and C iv emission lines. The lag profile in the blue wing of Ly α is consistent with virial motion, with longer lags dominating at lower velocities, and shorter lags at higher velocities. The C iv lag profile shows the signature of a thick rotating disk, with the shortest lags in the wings, local peaks at ±1500 km s −1 , and a local minimum at the line center. The other emission lines are dominated by broad absorption lines and blending with adjacent emission lines. These require detailed models, and will be presented in future work. 

 We present the first results from a 100-day Swift, NICER, and ground-based X-ray–UV–optical reverberation mapping campaign of the Narrow-line Seyfert 1 Mrk 335, when it was in an unprecedented low X-ray flux state. Despite dramatic suppression of the X-ray variability, we still observe UV–optical lags as expected from disk reverberation. Moreover, the UV–optical lags are consistent with archival observations when the X-ray luminosity was >10 times higher. Interestingly, both low- and high-flux states reveal UV–optical lags that are 6–11 times longer than expected from a thin disk. These long lags are often interpreted as due to contamination from the broad line region; however theu-band excess lag (containing the Balmer jump from the diffuse continuum) is less prevalent than in other active galactic nuclei. The Swift campaign showed a low X-ray-to-optical correlation (similar to previous campaigns), but NICER and ground-based monitoring continued for another 2 weeks, during which the optical rose to the highest level of the campaign, followed ∼10 days later by a sharp rise in X-rays. While the low X-ray countrate and relatively large systematic uncertainties in the NICER background make this measurement challenging, if the optical does lead X-rays in this flare, this indicates a departure from the zeroth-order reprocessing picture. If the optical flare is due to an increase in mass accretion rate, this occurs on much shorter than the viscous timescale. Alternatively, the optical could be responding to an intrinsic rise in X-rays that is initially hidden from our line of sight.

 

Abstract Quasar black hole masses are most commonly estimated using broad emission lines in single epoch spectra based on scaling relationships determined from reverberation mapping of small samples of low-redshift objects. Several effects have been identified requiring modifications to these scaling relationships, resulting in significant reductions of the black hole mass determinations at high redshift. Correcting these systematic biases is critical to understanding the relationships among black hole and host galaxy properties. We are completing a program using the Gemini North telescope, called the Gemini North Infrared Spectrograph (GNIRS) Distant Quasar Survey (DQS), that has produced rest-frame optical spectra of about 200 high-redshift quasars (z = 1.5–3.5). The GNIRS-DQS will produce new and improved ultraviolet-based black hole mass and accretion rate prescriptions, as well as new redshift prescriptions for velocity zero points of high-z quasars, necessary to measure feedback.